Bosonic stringlike behavior and the ultraviolet filtering of QCD 
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The gluonic action density is calculated in static mesons at finite temperature just below the 
deconfinement point. Our focus is to elucidate the role of vacuum ultraviolet fluctuations which 
are filtered using an improved smearing algorithm. In the intermediate source separation distance, 
where the free string picture poorly describes the flux tube width profile, we find upon reducing the 
vacuum action towards the classical instanton vacuum, the characteristics of the flux tube converge 
and compare favorably with the predictions of the free bosonic string. This result establishes a 
connection between the free string action and vacuum gauge fields and reveals the important role of 
ultraviolet physics in understanding the lattice data at this temperature scale. As a by-product of 
these calculations, we find the broadening of the QCD flux tube to be independent of the ultraviolet 
filtering at large distances. Our results exhibit a linearly divergent pattern in agreement with the 
string picture predictions. 

PACS numbers: 12.38.Gc, 12.38.Lg, 12.38.Aw 



I. INTRODUCTION 

In the dual superconductor scenario of quark confine- 
ment, the quantum chromodynamical (QCD) vacuum 
squeezes the color fields into a confining string dual to the 
Abrikosov line by the dual Mcissner effect. The string 
conjecture [l[ follows as an intuitive realization of this 
squeezed color field with the major objective of deriving 
the leading and subleading properties of the flux tube in 
the infrared region of confining gauge theories. This ef- 
fective description is expected to hold on distance scales 
larger than the intrinsic thickness of the flux tube 1 /T c Q 
in the rough phase of lattice gauge theories (LGT) . The 
linearly rising potential part arises from the classical con- 
figuration of the string, and the quantum fluctuations of 
the string lead to the presence of a long distance c/r 
term in the qq potential well known as the Liischer term. 
The existence of the subleading term has been verified in 
high precision measurements of Polyakov loop correlators 
in the SU(3) gauge group at zero temperature Q. The 
fluctuations of the string render an effective width for 
the flux-tube which grows logarithmically Q as the color 
sources are pulled apart. The logarithmic divergence has 
been verified in many lattice simulations corresponding 
to a variety of confining gauge models 0, [H-Q . 

At high temperature, higher-order gluonic modes are 
present. The corresponding free bosonic string predicts 
a new set of measurable thermal effects. These include 
a decrease in the effective string tension [8l-fl0j , a change 
in the pattern of the tube's growth in width from a log- 
arithmic divergence into a linear divergence (Tlj . and a 
non-constant width profile [ill, E3 along the qq line. 

Unlike the situation at zero temperature, the thermal 
behavior of the free string manifests only at source sepa- 



ration distance scales larger than what one expects nor- 
mally in the zero temperature regime [TTl - flij . The fact 
that the lattice data are poorly described by the free 
theory in the intermediate distance regime has been a 
subject of analytic and numerical studies which include 
higher-order terms of the effective string's action [l5|, [T(| 
into the corresponding partition function. The conse- 
quences of such an approach have been studied on the 
level of the qq potential [13, [H| and, recently, has been 
extended to the flux-tube width profile 0,01, Ell- Other 
studies investigate a possible finite intrinsic thickness of 



the QCD flux-tube [21 
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Apart from the linearly rising potential, the interest- 
ing physics of the effective confining string is mainly due 
to its quantum fluctuations. As we will see, remarkable 
features arise when the ultraviolet (UV) part of the fluc- 
tuations of the stringlike flux tube is filtered out for in- 
termediate quark separations at high temperatures. At 
this distance scale it is not yet clear if the deviations 
from the string picture are due to a non-Nambu-Goto 
action or the fact that a stringlike behavior has not yet 
set in. Thus, it is interesting to address this problem in 
a variant context by reporting an observation regarding 
the role played by the UV fluctuations of the vacuum in 
these discrepancies. We do this by tracking the response 
of the QCD vacuum, which is subject to UV filtering, to 
the presence of external static color sources. This work 
extends the region for which the free string picture is of 
utility. 

In the following we measure the gluonic action-density 
distribution by correlating an action density operator to 
Polyakov loop correlators. Measurements are taken on 
a set of SU(3) pure gauge configurations. The configu- 
rations are generated using the standard Wilson gauge 
action S w on two lattices of a spatial volume of 36 3 and 
temporal extents of N t = 10 and N t = 8, correspond- 
ing to temperatures T/T c w 0.8 and T/T c w 0.9 respec- 
tively. The simulations are performed for coupling value 
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(3 = 6.00. At this value the lattice spacing is a = 0.1 
fm to reproduce the standard value of the string tension 
sfB = 440 Mev The Monte Carlo updates are im- 

plemented with a pseudo-heatbath algorithm [23 1 using 
Fabricius-Haan and Kennedy- Pendelton (FHKP) 24l.[25| 
updating. Each update step consists of one heat bath and 
5 over-relaxations. The measurements are taken on 500 
bins separated with 2000 updating sweeps. Averaging 
inside each bin is performed by taking 5 measurements 
separated by 70 updating sweeps. This leads to a hier- 
archical integration, that is apart from updating the last 
time slice, similar to implementing a one-level Luscher 
Weisz (LW) algorithm 

The measurements arc taken after smoothing the 
gauge field by an over-improved stout-link smearing al- 
gorithm [26j . The value of the smearing parameters used 
are e = —0.25 and p^ = p = 0.06. Smoothing the gauge 
field reduces the action towards the action minimum or 
the classical instanton solution (SJ. The UV character- 
istics of the gauge fields can be characterized in terms 
of Dirac eigcnmodcs. For example, the number of over- 
improved stout-link smearing sweeps used here has been 
calibrated to a given spectral cut-off A cu toff in the spec- 
tral representation of the Dirac operator [28|. The mea- 
surements are taken on sets of smeared gauge configura- 
tions with increasing levels of smearing. This way we are 
able to set the limit where the QCD vacuum response to 
the presence of an external static color sources asymp- 
totically approaches the low energy free effective theory 
behavior. 



II. QUARK-ANTIQUARK POTENTIAL 

At fixed temperature T, the Monte Carlo evaluation 
of the quark-antiquark potential at each R is calculated 
through the Polyakov loop correlators 



A(Ci,C2 = R) — 0, and periodic boundary condition 
along the time direction X((i = 0, £2) = -X"(Ci = Lt, (2), 
Lt = if. S is the string action in the physical gauge [3| 



s[x] = a- + - j o da y o d(, 



(2.4) 



The action decomposes into the classical configuration, 
the fluctuation part, and the string higher-order self in- 
teractions. A leading order approximation can be made 
by neglecting the self-interaction terms. £ = (£1,^2) is 
a complex paramctrization of the world sheet, such that 
G [-R/2,R/2],& G [-L T /2,L T /2], and r = ^ is 
the modular parameter of the cylinder. Solving the path 
integral of Eq. ([2~4]l [29j . and using Eq. ([2TT]) , the quark 
anti-quark potential reads 



V{R,T) = 2T\ogr ] 



2TR 

rj is the Dcdckind eta function 



(jR + p{T). (2.5) 



t?(t) 



(2.6) 



and /i(T) is a renormalization parameter. 

The numerical evaluation of the quark anti-quark po- 
tential, Eq. (|2.1|) . using a four dimensional smearing 
scheme leads to a systematic ambiguity in regard to the 
transfer matrix interpretation which allows one to iden- 
tify the expectation values of the Polyakov loop corre- 
lators with exp(— V(R, T)/T). We recourse, instead, to 
three dimensional smearing keeping the temporal-links 
unsmeared. The same smearing parameters as above are 
used. 



V 2Q = J d[U]P(0)P\R)eM-S w ), 
= eM-V(R,T)/T). 
with the Polyakov loop given by 



(2.1) 



P(n) = ^Tr 



»if=i 



(2-2) 



In the string picture, the Polyakov loop correlator as- 
sumes the functional form of the partition function of the 
two-dimensional bosonic string 



(P(0)pt(iJ))= f[DX}exp(-S(X)). 
Jc 

The vector ^(£1^2) maps the region C C R 2 
R 4 , with Dirichlct boundary condition A (£1,(2 = 



(2-3) 
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FIG. 1: The temporal link Ui is updated based on the 
neighboring links. The shaded area represents the 3D spatial 
smeared lattice. The heat bath starts from smooth spatial 
links. 
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(a) 3D smearing. 
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(b) 4D smearing. 



FIG. 2: The quark-antiquark potential measured at each depicted smearing level for three and four-dimensional smearing, 
the lines correspond to the string picture predictions of Eq. (|2.5p . The standard value of the string tension is used. 



Our approach is as follows: 

• We start with 500 configurations. Each configura- 
tion is smeared only altering the spatial directions 
of the links, using spatially-oriented staples. This 
is depicted by shading of spatial links in Fig. [TJ 

• For every 3D smeared configuration, an update 
sweep is applied. Each update sweep consists of 
a heat-bath step and four overrelaxation steps on 
the lattice. Four or more overrelaxation steps pro- 
vide the same result within errors. Each step pro- 
ceeds as an update of every single link based on its 
neighbors, effectively one at a time as one sweeps 
across the lattice in all four directions; spatial and 
temporal as indicated in Fig. [TJ 

• An update sweep and subsequent measurement is 
repeated three more times on each configuration of 
the ensemble. This results in four measurements 
one over each of the four ensembles, where the en- 
sembles differ in how many update sweeps they 
have had after spatial smearing, i.e, 1 to 4. 

The measurements proceed as follows: 

• The temporal links in each of theses newly created 
configurations arc integrated out using a source Q 
sum of staples. The temporal link variables Ut are 
replaced with the new link variable 



using the numerical link-integration method of 
Ref. L30i. It should be noted, that unlike the 



fdUUe- Tr ^ ui + u Q f ) 
/d[/e- Tr W f,t + t "3 t ) ' 



(2.7) 



link integration implementation of Ref. [31[ where 
pseudo-heat bath hits are performed only on the 
temporal links, the pseudo-heat bath updates we 
described above are performed on both the smeared 
space oriented links and the time oriented links. 
The temporal links, however, are integrated out by 
the numerical evaluation of the equivalent contour 
integral of Eq. (|2.7[) as detailed in Ref. [3(J ■ 

• Finally, the Polyakov loop correlators are calcu- 
lated on each of the four one update sweep sep- 
arated configurations. The result is averaged and 
binned as a single jackknifc entry to avoid artificial 
error reduction. The 500 decorrelated bins are then 
averaged. 

The Monte Carlo update step starts from a low action 
configuration in the spatial directions due to smearing. 
The above described update procedure brings in a newly 
updated time-link such that the effects of local action re- 
duction that was only in the spatial torus takes place in 
the four-dimensional lattice. In this way, the UV filtering 
is implemented keeping the integration over the path in- 
tegral Eq. (|2.1[) systematic, thus, preserving the transfer 
matrix interpretation. 

In Fig. [5] the value of the potential measured on various 
levels of spatially smeared configurations, normalized to 
its value at R = 1.2 fin, are plotted. Fig. 2(a) shows the 
numerical behavior of the data using the above described 
3D smeared heat bath/overrelaxation driven updates. 
On the other hand, Fig. 2(b) shows the corresponding 
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numerical behavior of the data measured on standard 
four-dimensional smeared configurations. The data cor- 
responding to the unsmeared lattice and the string model 
predictions of Eq. (|2.5p at T/T c ps 0.9 are also included. 

The discrepancies between the unsmeared lattice data 
and the free string model occur in the intermediate dis- 
tances R < 1 fm. The numerical results for the quark- 
antiquark potential evaluated on the 3D smeared up- 
dated configurations show an interesting behavior with 
respect to the number of smearing sweeps (see Fig. 2(a)). 
The data at large distances show no response to the fil- 
tering of the UV fluctuations of the gauge field. The 
data at intermediate separation distances converge for a 
large number of smearing sweeps. Moreover, the data 
approach the free string model predictions. 

It is interesting to compare these results to those ob- 
tained from the four-dimensional smearing illustrated in 
Fig. 2(b), where 20 and 100 sweeps of smearing are com- 
pared. The results for 100 sweeps of smearing coincide 
very well with the convergence toward the string model 
predictions observed in Fig. 2(a). Thus, the AD smearing 
approach can be used as an efficient method for explor- 
ing the more demanding three-point functions required 
to determine the distribution of gluon flux, and this is 
used in the following. Similar results are observed at 
T/T c ps 0.8. 



III. THE GLUONIC PROFILE 

The transverse degrees of freedom of the string-like 
flux tube render an effective width for the tube. The 
mean-square width of the free bosonic string is defined 
as 

W 2 (£;t)^ (A 2 (£;t)) 

J c [DX]X 2 eM~S[X}) 
J c [DX]eM~S[X}) ■ [ ■ ' 

The above integral can be solved analytically [Tl], . 
The width of the tube in D dimensions reads 



^ 2 (5,t) 



D 



27TCT 



■log 



R 



Ro(0 



D 



27TCT 



log 



e 2 {n£/R;T) 



(3.2) 

where 8 are Jacobi elliptic functions, and i?o(£) is the 
UV cutoff which has been generalized to be dependent on 
distances from the sources. This solution gives the mean 
square width at all the planes transverse to the quark- 
antiquark line, and hence, is describing the topological 
shape of the fluctuating flux tube and its dependence on 
the temperature as well as its evolution with the increase 
of color source separation. In the above formula, we as- 
sume a dependence of the ultraviolet cutoff, i?o, on the 
position of the transverse planes, since as we will see, this 
quantity assumes different values near the quark sources 
as we fit the above formula to lattice data. 



Using a modular transform r — > — 1/r (Til [lj|, 
Eq. p.2[) , at R » i, in four dimensions (2(| becomes 



^(R/2;T) = ±lo g ( 2 -^ 
ira \ 1 



±-TR-—e- 2 * RT 
2a na 



(3.3) 

which indicates linear growth of the tube's width at large 
distance. 
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FIG. 3: Plot of the density distribution C(z = R/2, xo,y) in 
the quark plane, xo, at the center of the tube, R/2, for source 
separation R = 0.9 fm, T/T c PS 0.9. 



The width of the action density of the free bosonic 
string can be compared to the width of the action den- 
sity of the corresponding flux tube of the lattice gauge 
theory. After constructing the color-averaged infinitely- 
heavy static-mesonic state 

subsequent measurement by an action density operator 
i (E 2 ~B 2 ) is taken at each point of the three-dimensional 
torus at each corresponding Euclidean time slice for ev- 
ery source configuration. The action density operator is 
constructed via a highly- improved Q(a ) three- loop im- 
proved lattice field-strength tensor |32j. The measure- 
ments taken are averaged over the time slices. A scalar 
field that characterizes the gluonic action-density distri- 
bution field can be then measured using the definition (33j 



C(p;ri,r 2 ) 



(V2 Q (rur 2 )S(p)) 
'V 2Q (r u r 2 )) (S(p)Y 



(3.4) 



where < > denotes averaging over configurations 

and lattice symmetries, and the vector p refers to the 
spatial position of the flux probe with respect to some 
origin. To further suppress the statistical fluctuations, 
the density distributions have been symmetrised around 
all the symmetry planes of the tube. 
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Unless otherwise indicated, the measurements pre- 
sented through out this section are taken on the 4D 
smeared configurations. The decrease mC(y) with the in- 
crease of the smearing sweeps is depicted in Fig. [3] The 
action-density asymptotically converges to a minimum 
for values around n sw = 80 to 100 sweeps of smearing. 
Cluster decomposition of the operators leads to C — > 1 
away from the quarks. 

A measurement of the width of the flux-tube's action 
density may be taken through fitting the density distri- 
bution C(p(z, r, 9)), Eq. (|3.4j) . to a Gaussian of the form 



C(z; r, 6) = 1 - a exp[-r 2 /^ 2 (^)] 



(3.5) 



with r 2 = x 2 + y 2 in each selected transverse plane 
p(z; r, 9) to quark axis z by making use of the cylindrical 
symmetry of the tube. 

The mean square width of the flux tube is defined as 
the second moment of the flux density with respect to 
the central line connecting the two quarks, z, 



L0 2 (Z) 



J drr 3 



7" 2 (*)] 



J drr exp[— r 2 /ui 2 (z)] 



(3.6) 



with z = or z = R denoting the position of the quark 
source. The width is normalized with respect to the ul- 
traviolet cutoff i?o (0 of Eq. (|3.2[) according to 

oj 2 n (z)=oj 2 (z) + —log(Ro(z)). (3.7) 
7r a 

The measured values of the mean-square width in the 
middle plane z = R/2 of the tube versus the source 
separation arc plotted in Fig. @] Similar plots at three 
consecutive transverse planes z = 2, z = 3, and z = 4 
to the line joining the two color sources are illustrated 
in Figs. [5] and [6] i?o(£) has been measured for each 
smearing level and £ value by fitting Eq. (|3.2p to data 
points having R > 1 fm. Good \ 2 is obtained. 

At large distances, in the middle plane of the flux tube, 
the tube shows a width broadening pattern for increas- 
ing R that docs not depend strongly on the corresponding 
smearing level. The data at large distances are increas- 
ing linearly in agreement with the string model predic- 
tions @, O, H(j. The UV effects, on the other hand, are 
manifest in the data points at shorter distances. The 
width of the flux-tube measured on the lowest smearing 
level, where the short distance physics is best preserved, 
is poorly described by the free bosonic string model at 
short distances. As higher smearing levels are considered, 
the subsequent removal of the short distance physics from 
the gauge sector regulates the fast rate of growth of the 
flux tube width. However, this does not continue uncon- 
trollably. The data ultimately converge near 100 sweeps 
of smearing, in accord with the saturation in the action 
density of Fig. [3] and the potential in Fig. 2(a). More- 
over, the UV- filtered results converge to the free string 
predictions. 




1 1.1 

R(fm) 



FIG. 4: The mean square width of the flux tube ui\{z = R/2) 
in the middle plane between the quarks. The lattice data, 
corresponding to the action density minimization, approach 
the string model predictions at short distances. At large dis- 
tances the predicted linear divergence of the flux tube width 
is manifest in lattice data. 



Table Q] summarizes the measured Xdof f° r ^ s of 
Eq. (|3.2I) for the fit range R > 0.5 fm. With the increase 
of the number of smearing sweeps the returned values of 
the Xdof improves and becomes stable near the regime of 
the action saturation of approximately 100 sweeps. 



TABLE I: The returned Xdof by the fit of lattice data for mea- 
surement on the flux tube width in the middle plane between 
the quark anti-quark to the effective string model predictions 
of Eq. (|3.2[) . The lattice data correspond to smearing levels 
from n sw = 40 to n sm = 120 



No.Sweeps 40 60 80 100 120 



Xdof 



3.2 1.6 1.20 0.98 0.96 



To clarify this point further, we investigate the re- 
sponse of the QCD vacuum to the presence of infinitely 
heavy sources that are not constructed using smeared 
temporal links. Instead, the Polyakov lines are evaluated 
in the calculations of the flux strength, Eq. (|3.4p . using 
the numerical link integration [3lj procedure of Ref. [30j ] 
for noise reduction. The temporal links have not been 
smeared in the evaluation of Polyakov loops, rather the 
temporal links have been integrated out. This time, the 
Polyakov loops are taken from the unsmcarcd configura- 
tions and correlated with the smeared action density. For 
our analysis, performed on 500 configurations, we observe 
that the data corresponding to the width profile of the 
flux tube measured on high levels of vacuum UV filtering 
do display similar behavior to the results in Fig.[5J Fig. [5] 
presents results for the plane z = 3. Again the results 
systematically converge with a large number of smearing 
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FIG. 6: Same as Fig[5]for LJn(z) measured at the plane z = 3. 
The Polyakov lines are evaluated after integrating out the 
time links. The legend indicates the number of smearing 
sweeps applied before action measurements are taken. 
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FIG. 5: The normalized width of the flux tube to^(z) versus 
qq separations measured in the planes (a) z = 2, (b) z — 
3, (c) z = 4, at T/T c « 0.9. The coordinates z are lattice 
coordinates (lattice units) and are measured from the quark 
position z — 0. The line denotes the one parameter string 
model, Eq. (|3.2[) . fit to lattice data for R > 1 fm. The numbers 
in the legend denote the number of smearing sweeps. 



sweeps and approach the string model predictions. Note 
that the evaluation of the correlation function using this 
method involves a three point correlation function which 
becomes noisy at large distances. 

The shape aspects of the fluctuating free string are 
contained mainly in the second term of Eq. (|3.2|) and 
can be isolated by considering the difference in the mean 
square width at a given plane with respect to the central 
plane, Sui 2 = w 2 (z) — uj 2 (zq). The measured value of Rq, 
however, depends on the corresponding plane at which 
the lattice data is fit to Eq. (|3 . 2[) . The value of the fit 
parameter Rq is fixed for each plane using lattice data at 
large separations R = 1.1 fm and R = 1.2 fm. 

The changes in Rq with respect to central plane is plot- 
ted separately in Fig. [7] The measurements performed 
at adjacent planes of the flux-tube reveal that this quan- 
tity varies along different planes orthogonal to the string. 
Indeed, we don't get a perfect match with the free-string 
profile unless we take such changes into account. The 
increase in the value of the UV cutoff Rq is mostly ob- 
vious near the quarks and may indicate the importance 
of string self interactions near the boundary or the string 
interaction with the quark source itself. These interac- 
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FIG. 7: The measured change in the ultraviolet cutoff Ro 
along the flux tube. The smearing effect is small relative to 
the nontrivial dependence of Ro on z. 
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of generalized string model predictions with the corre- 
sponding change in the flux tube's width on the lattice 
(obtained from fits of Eq. (|3.4[l ) is shown in Fig. [5] for a 
source separation in the intermediate distance R = 0.8 
fm, R = 0.9 fm and at large distance R = 1.2 fm. 

Lattice data for each gauge smoothing level is also 
depicted in Fig. [5] The jackknifc uncertainties associ- 
ated with the change in the mean-squared width of the 
tube reveal correlated errors between the adjacent planes. 
Only subtle changes are observed in the tube's width 
along the transverse planes with respect to the central 
plane for the analysis performed on smeared gauge con- 
figurations of the lowest smearing level n sw = 40 sweeps. 
The tube tends to exhibit larger curvatures as higher- 
levels of gauge smoothing arc considered. At large val- 
ues of UV filtering, the tube profile converges and ap- 
proaches the geometrical shape of the free-bosonic string. 
At large distances, on the other hand, the flux tube dis- 
plays a curved width profile which compares well with 
the bosonic string profile and is not affected by smear- 
ing, as is evident in Fig. [8] at R = 1.2 fm. The string 
model curve has been calculated based on the changes in 
the second term of Eq. 
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FIG. 8: The change of the tube's width Suj 2 = uj 2 (z) - uj 2 (z ) 
measured from the central plane for the depicted qq separa- 
tions. The smearing level of the lattice data is illustrated. 
The line denotes the width difference Suj 2 according to the 
string model Eq. ||32J|. P = 6, T/T c w 0.9. 



tions are switched-off in the Nambu-Goto (NG) action 
and would appear when fitting the lattice data to the 
free string profile as a variation in the value of the UV 
cutoff along the tube. 

This variation in R has been included as part of the 
generalized string model solution of Eq. (|3.2j) illustrated 
by the curves in the following figures. A comparison 



The presence of a pair of static external sources in 
the QCD vacuum induces a response of an effective free 
bosonic string for source separations in the intermedi- 
ate separation region, provided the short distance vac- 
uum fluctuations are filtered out. The flux tube, mea- 
sured as a correlation between the mcsonic operator and 
the vacuum action density, is found to exhibit a broad- 
ening pattern and a transverse structure similar to the 
free-bosonic string for measurements taken near the sat- 
uration in action density minimization under smearing 
even at intermediate distances. At large distances, the 
UV fluctuations do not affect the tube growth, which 
exhibits a linear divergent pattern consistent with the 
string model predictions. 

In carrying out the qq potential calculations, we intro- 
duced a novel method for studying the effects of UV fil- 
tering of the QCD vacuum. This method avoids the am- 
biguities of performing unsystematic integrations due to 
smearing the temporal links, thus, preserving the transfer 
matrix interpretation. Instead three-dimensional spatial 
smearing is combined with single pseudo-heatbath driven 
updates. The numerical data of the qq potential obtained 
this way converge towards the string model predictions 
at large number of smearing sweeps. A comparison with 
four-dimensional smearing results reveals that any sys- 
tematic effects associated with smearing the temporal 
links are subtle. 

The analysis performed at short distances provides an 
extension of the QCD vacua where the free string pic- 
ture is of utility. The infrared region of the UV filtered 
vacuum can be described merely on the basis of a free 
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string picture in the intermediate distances as well as 
large distances. This fact is relevant and complements 
recent investigations including higher-order self interac- 
tions to match lattice results. 
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